Deformation of Small Compressed Droplets 
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We investigate the elastic properties of small droplets under compression. The compression of 
a bubble by two parallel plates is solved exactly and it is shown that a lowest-order expansion of 
the solution reduces to a form similar to that obtained by Morse and Witten. Other systems are 
studied numerically and results for configurations involving between 2 and 20 compressing planes 
are presented. It is found that the response to compression depends on the number of planes. The 
shear modulus is also calculated for common lattices and the stability crossover between f.c.c. and 
b.c.c. is discussed. 



PACS numbers: 82.70.Kj, 81.40. Jj, 62.20.Dc 
I. INTRODUCTION 

Emulsions consist of a mixture of two immiscible flu- 
ids, one of which, generally an oil, is dispersed as small 
droplets in the continuous phase of the other fluid, gener- 
ally water. The interfaces are stabilized by a surfactant, 
preventing coalescence. Emulsions are materials with 
quite unusual properties: despite being comprised solely 
of fluids, they become elastic solids when the droplets 
are compressed to a large enough volume fraction tp by 
extracting the continuous phase by the application of an 
osmotic pressure II. The origin of the elasticity is the in- 
tcrfacial energy of the droplets. At low volume fractions, 
the surface tension a ensures that the droplets are spher- 
ical in shape. However, at higher ip, the packing con- 
straints force the droplets to deform, thus storing energy. 
Consequently, it is a fundamental issue to determine the 
increase of the surface area of a droplet resulting from an 
arbitrary deformation. 

It is generally believed that the application of a shear 
strain to a compressed emulsion causes the droplets to 
further deform, thereby storing more energy [|]J|^,^,||. 
According to this picture, the scale of both the osmotic 
pressure, and the elastic shear modulus G is set by the 
Laplace pressure of the droplets, (2a /R), where R is the 
radius of an undeformed droplet. Recent experiments j(| 
on the elastic properties of compressed monodisperse 
emulsions of silicone-oil in water confirm the role of the 
Laplace pressure: the experimental values of G and LT for 
different droplet sizes, when scaled by (cr/R), collapse on 
a single curve |J . 

Compression modifies the elastic properties of a dis- 
ordered emulsion, changing the response from liquid-like 
to solid-like. The onset of solid-like behavior is gradual; 
both II and G increase smoothly from zero as the system 
is compressed above tp c 0.64 the volume fraction 
at which disordered monodisperse droplets are first de- 
formed 0. For polydisperse emulsions, ip c is larger ||, 
reflecting more efficient packing. Two separate elastic re- 
sponses of the droplets are therefore of interest: to direct 
compression and to shear. 



The increase of the surface area of a droplet under com- 
pression comes from the relatively low cost of deforma- 
tion compared to the compressibility of the internal fluid. 
Thus, it is very reasonable to assume that the droplets 
are composed of an incompressible fluid, and may be con- 
sequently considered to have a fixed volume. The effect 
of gravity may be nullified by density matching of the oil 
and water, but in any case would be tiny because of the 
large surface tension and small droplet size (typically of 
a /itm) of the experimental systems. 

When the droplets are compressed just beyond the 
point they start to deform, i.e. (p ~> (p c , they form small 
facets separated by thin surfactant-water-surfactant 
films. Each of the facets created on the droplet surface 
has the effect of locally decreasing the droplet surface, 
but volume conservation makes the droplet globally in- 
crease its surface, and thus its energy. Due to such non- 
local effects, the determination of the droplet response 
to compression requires that it be considered as a whole. 
Consequently, a complete analytic solution is limited to a 
few simple cases, one of which is a drop squeezed between 
two parallel plates, to be presented in Section [II. The 



response to shear is still more difficult, but fortunately, 
numerical tools capable of handling such problems have 
recently become available || . 

Upon more compression, droplets gradually take the 
shape of a rounded polyhedron. At this point, the water 
is contained within the connected network of voids left 
between the oil droplets. This network is made of thin 
veins, along which the facets meet, and which connect 
larger regions located at the corners of the polyhedra. As 
ip — ► 1, i.e. in the so-called dry foam limit, mechanical 
equilibrium imposes constraints on these regions which 
were first discussed by Plateau 10 1. In particular, ex- 
actly three films meet at equal angles (120°) along the 
veins which become the edges of the polyhedra. Addi- 
tionally, the corner regions reduce to a point at which 
four edges meet at equal tetrahedral angles. These rules 
do not apply until most of the continuous phase has been 
extracted. 

The elasticity of emulsions has been the object of nu- 
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merous investigations over the last decade. The first the- 
oretical attempts were for ordered lattices in two 
dimensions, which we shall briefly treat in Section ||, 
mainly in order to contrast them with three-dimensional 
systems. It should be noted that two-dimensional sys- 
tems are important in their own right, both theoretically 
and experimentally. For example, the effects of polydis- 
persity and disorder are largely studied in two dimen- 
sions 

For real (three-dimensional) bubbles, the only analyti- 
cal expression [Q for the behavior of a droplet surface un- 
der compression involves the application of an infinites- 
imal force at a single point. A simplified approach con- 
sists in assuming that each droplet becomes a truncated 
sphere under compression Using each of these ap- 
proaches, the shear modulus of a simple cubic (s.c.) lat- 
tice has been derived f|. These studies predict that the 
onset of the shear modulus at ip c is discontinuous, in 
contrast with the smooth quasi-linear increase of G at Lp c 
found experimentally for disordered emulsions. 

The relative importance of the disorder of the emul- 
sions, its polydispersity, and the response of the droplets 
themselves is the object of current research 14 1. This pa- 
per is concerned with the last: we present the results of 
an investigation of the response of an individual droplet 
to deformation. In Section [n we present results for the 2- 
dimensional case. Section III derives exact results for the 
compression of a droplet by two parallel plates. In the 
limit of infinitesimal compression, our results are com- 
pared with Morse and Witten's expression Q . For a large 
range of compressions, we show that the energy has an 
anharmonic power-law response to compression. In Sec- 
tion IV, we present numerical results for the compression 



of a droplet in several local environments using Brakke's 
Surface Evolver (SE) software || . We again demonstrate 
that the energy of the squeezed droplet is anharmonic, 
with the power depending on the number of facets. Sec- 
tion ^ describes the results obtained for a selection of 
space-filling structures at the dry foam limit. Section VI 
contains numerical calculations of the droplet response to 
a shear deformation. We discuss our results and conclude 



in Section VII 



II. RESPONSE IN TWO DIMENSIONS 

Two-dimensional systems have characteristics not 
found in three-dimensions. For example, monodisperse 
arrays of disks order easily Jig ] in an hexagonal array 
of surface fraction ip c — tt/(2\/3) « 0.9069. The elastic 
properties of two-dimensional ordered systems have been 
discussed extensively by Princen |}],|l6). The response to 
deformation is also simpler in two dimensions. A min- 
imum free surface (a free curve in the present case) is 
characterized by a uniform pressure or, equivalently by 
a constant mean curvature. In two dimensions, the sur- 
face is parameterized by only one radius of curvature and 
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FIG. 1. The excess energy per facet of a two-dimensional 
droplet as a function of compression. The curves represent, 
from right to left, the tessellating cases n = 3, 4, and 6. The 
insert shows the same data on a logarithmic scale. 



the minimum free surface is therefore always an arc of a 
circle. 

Let us consider the deformation of a two-dimensional 
droplet of constant area confined inside a regular polygon 
of n sides. Only for n — 3 ,4, and 6 can such polygons 
be tessellated, but we shall consider general n, in partic- 
ular to contrast the two- and three-dimensional results. 
The ratio of droplet area to polygon area is denoted by if. 
Only for the tessellating cases will ip represent the surface 
covering fraction of the corresponding lattice. In the fol- 
lowing, we assume that the droplet is non-wetting so that 
the contact angles are zero. At (p c = — cot — , the circular 
droplet is undeformed and touches the polygon at exactly 
n points, being the midpoints of the polygon's edges. If 
the polygon is now shrunk uniformly, the droplet will dis- 
tort, and its new minimum-perimeter shape will consist 
of small facets joined by circular arcs. As if increases 
above (/J c , the length of the flat portions increases and 
the radius of curvature of the arcs decreases, reflecting 
increased droplet pressure. At <p = 1, the radius of cur- 
vature of the arcs becomes zero (infinite pressure), and 
the bubble takes on the shape of the polygonal cell. 

The interfacial energy of the droplets is the product 
of their line tension a and the perimeter length £. It is 
convenient to measure the excess energy by the following 
dimensionless quantity 



2ttR 
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(1) 



where R is the radius of the undeformed circular droplet. 
Similarly, the degree of compression will be measured by 
the following dimensionless displacement 
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R-h 
R ' 



(2) 



where h is the perpendicular distance from the facets 
to the center of the droplet. £ is simply related to the 
volume fraction through £ = 1 — y/tp c /tp. For any regular 
polygonal cell, e can be shown to be, for ip > ip c , 



(3) 



For 93 slightly greater than ip c , we may expand the pre- 
vious expression in terms of Up — ip c ), and get 
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2(1 - ifc) 



(4) 



Therefore, for small compression, the energy is harmonic, 
and the prefactor depends on the number of faces of the 
polygon through <p c . This is shown in Fig. |l| where the 
excess energy per facet is plotted for the tessellating cases 
n = 3, 4 and 6. The end of each curve corresponds to 
ip — ► 1 for each lattice. The logarithmic insert shows the 
wide range over which harmonicity holds as well as the 
n dependence of the prefactor expressed by the different 
intercepts. 

In two dimensions, the scaled osmotic pressure can be 
obtained from 



U/(a/R) = 2ip' 



de 
dip' 



and we find, using Eq. 0, 



U/(a/R) = ( 2- 

\9c 



1 - <Pc 

l-ip 



(5) 
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is 



This last expression, which depends on n through ip, 
valid for all tessellations in two dimensions. For (p > 
II can be shown to increase linearly with Up — tp c ). 

In the model proposed by Princen Jl],[l6|], droplets are 
monodisperse cylindrical objects packed in an hexagonal 
array (n = 6). In the dry foam limit, this configuration 
satisfies Plateau's rules for packing in two dimensions. 
For ip > ip c , the osmotic pressure is given by Eq. ^ with 
n = 6 Jliyiq] . For the same range, the static shear mod- 
ulus was found to obey Q 



G/(a/R) = 0.525^. 



(7) 



Thus, this model predicts a discontinuity of G at tp c , its 
value jumping from zero to 0.525^/^. The onset of the 
shear modulus is thus very sharp at ip c , and we shall 
demonstrate that this elastic response is intimately as- 
sociated with the response of the droplet potential at 




FIG. 2. A droplet of initial radius R being deformed by 2 
parallel plates located at z = —h and z = +h. The equator is 
located at z = and has a radius r e . Both facets at z — ±h 
are circular and have a radius r/. The droplet shape is a cubic 
spline drawn for visualization purposes only. 



III. COMPRESSION BY TWO PLATES 

The evaluation of the excess energy stored in an arbi- 
trary surface deformation of a compressed droplet is a dif- 
ficult problem. A simple three-dimensional case consists 
of a droplet of radius R compressed between two parallel 
planes, each located at a distance h from the center of the 
droplet (cf. Fig.||). We shall measure compression by the 
dimensionless ratio £ defined as above. At small compres- 
sion d£, small circular facets of radius ry appear where the 
droplet touches each plate. Naively, the resulting force 
dF on each facet can be estimated by assuming that the 
radius of the droplet, and hence the Laplace pressure, 
remains unchanged. Therefore, dF w (2a / R)dS , where 
dS is the surface of the flattened facet. To lowest order 
in the deformation, 



dS PS 27ri? 2 d£, 



(8) 



so that dF « AiraRd^. Thus, to lowest order, the re- 
sponse of a droplet compressed between 2 parallel plates 
is found to be identical to the compression of a repulsive 
harmonic spring of spring constant Ana. However, this 
simple derivation leads to a wrong answer, as we now 
show. 

Because of its azimuthal symmetry, we may express 
this problem as a one-dimensional problem, that of find- 
ing the solid of revolution of a constant volume with min- 
imum surface area. The free surface of this solid is given 
by a curve r[z) rotated about the z-axis. The mathe- 
matical aspects of the present problem, such as stability 
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and the existence of a solution as a function of contact 
angles, have recently been discussed in detail by various 
authors @||Jl9g]. 

Using the Euler-Lagrange formalism, we minimize the 
total droplet surface, A, which is given by 



A = 2nrf + 2 1 + r\ dz 



(9) 



where r z = dr/dz. The minimization is done under the 
constraint of constant volume: 
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(10) 



which is introduced through the use of a Lagrange mul- 
tiplier A in a function C, 



£ = 
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2r^T 



- Ar 



(11) 



to be minimized. Since C = C(r,r z ) is independent of 
z, one can use the integrated form of the Euler-Lagrange 
equation 21 



dC 

dr z 



c = c, 



(12) 



where C is an integration constant. We re-express the 
unknown constants A and C using the following boundary 
conditions: 



r z \ rf =oo, 
rzl =0, 



(13a) 
(13b) 



where r e is the maximum value of r located at the equator 
(cf. Fig. ||). The first condition sets the contact angle 
at the facet, which should be zero if no long-range at- 
tractive forces are present |j],[l3| . The second comes from 
symmetry and macroscopic smoothness at the equator. 

It will be convenient to use the following dimensionlcss 
variables 



P = r/r e , 
Pf = r fl r e, 
C = z/r e , 

and, including these changes, Eq. fL2| reads 
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— P(l-P) 
p\ ~ ~ 



P - Pf 

Solving for z and integrating, we get 
z = —h + r e I(p), 



(14a) 
(14b) 
(14c) 



(15) 



(16) 



where p goes from pj to 1, and I(p) is an integral that 
can be solved at least numerically, and which is given by 
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FIG. 3. Shape of a droplet compressed between two paral- 
lel plates for different plate displacements: £ = 0.000, 0.035, 
0.107, 0.199, 0.302, 0.411, 0.526, 0.647. 
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PS ^J(l — x 2 )(x 2 — pj) 



dx. 



(17) 



Setting z = in Eq. [16| we obtain a relation between h/r e 
and pf: 

h = r e h (18) 

where I\ = 1(1)- This allows us to re-express Eq. [l6| as 

z = -h [1 - I( P )/h] . (19) 

We note here that both I(p) and I\ are functions of 
Pf, which in turn is related to h through the volume con- 
straint, Eq. |l^. Using Eqs. |l^ and [l8| h can be related 
to pf through the following non-trivial relation 
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3Ji 
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where J\ is another integral defined by 

, f 1 x\x^p]) 
•h — 



ht j{l-x*){x*-pj) 



dx. 



(20) 



(21) 



At this point the problem is fully determined. Given 
the size of the facets ry , the plate separation 2h is deter- 
mined by Eq. ^O], and the shape of the droplet by Eq. [nl 
The total droplet surface is obtained by combining Eqs]^, 
Il5|, and 18, to get 
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(22) 
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FIG. 4. Relative excess energy per facet for a single droplet 
as a function of the dimensionless displacement £. Curves 
are, from top to bottom: the truncated sphere approxima- 
tion (dot-dashed), the exact solution (solid), and the poten- 
tial predicted by Eq. ^ (dotted). Data points (o) are results 
obtained from SE. The absolute difference A in e between the 
semi-circular surface approximation and the exact solution is 
shown as an insert. 



where K\ is yet another integral defined by 



ps ^{1 - x*)(x* - pj) 



dx. 



(23) 



The integrals I± , J\ , and K\ are elliptic integrals which 
we can solve numerically, after removing the divergence 
at x = 1 by direct analytical integration. The solution 
of Eq. |l9], for various values of h, leads to the different 
droplet shapes shown in Fig. |[ 

It is convenient to define the (dimensionless) relative 
droplet excess surface energy by 



A/(4nR 2 



1. 



(24) 



Figure ^ shows the relative excess energy as derived from 
Eq. ^2| It is interesting to compare our solution to differ- 
ent estimates, which model the deformed droplets in dif- 
ferent ways. For example, we consider the relative excess 
energy of a body of equal volume consisting of a cylinder 
of radius surrounded by a semi-circle of revolution of 
radius h. This solution has a zero contact angle for all 
compressions and becomes exact in the two displacement 
limits £ — > 0, and £ — > 1. Given the scale of Fig. ||, it 
is not possible to distinguish the results of this estimate 
from the exact solution, and we thus show the difference 
as an insert. A second estimate which has been used in 
the case of compression by multiple planes |^Jl^| consists 



in considering a truncated sphere having the same vol- 
ume. Results from this estimate are also shown in Fig. |J. 
The force required to compress the plates is a better ac- 
curacy indicator, and forces derived from each solution 
will be compared below. 

The integrals 1\ , J\ and K\ can be solved analytically 
when pf = 0, i.e., for the undcformed sphere. Let us 
consider small compressions, and expand the integrals 
about pf — 0. For K\, for instance, we obtain 



(25) 




which becomes, after performing the integrals, 



? 4- r 

2 2 



ln(2) -I- \n{ Pf ) 



(26) 



In this limit one finds, after expanding Eq.^2| in a similar 
way, 



IX* ( \n{X) + ~ 



where 



X = (p f /2f 



(27) 



(28) 



has been introduced for simplicity. Using similar expan- 
sions in Eq. pel, one obtains 



£ ps -2X\n(X). 



(29) 



Due to this non-trivial relation, we cannot express e in 
terms of £. It is possible, however, to combine Eqs. |27] 
and [2^ to obtain an expression for the dimensionless 
force / 



de 



AX. 



(30) 



Combining the last equations, we may re-express the 
excess energy in terms of /: 



Hf)] 



(31) 



where the constant k = 2 ln(2) — |. Morse and Witten jij] 
obtained a similar result by considering, to lowest-order, 
a point perturbation of a sphere. In order to test the 
range of validity of our expansion, we solve Eq. |3l] nu- 
merically to obtain e(£). This is done by transforming 
Eq. |3l] into a self-consistent integral equation 



f(0 



4j[ttx)d^ 
k ~ m(/(£)) 



1/2 



(32) 



from which e(£) = J Q f(x)dx is obtained after a few tens 
of iterations only The results are shown by the dotted 
lines on Figures Q and [| 
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FIG. 5. A comparison of the £ dependence of the dimen- 
sionless force ratio //£ for different solutions. Curves are, top 
to bottom: the truncated sphere approximation (dot-dashed), 
the semi-circular surface approximation (dashed), the exact 
solution (solid), and its lowest-order expansion (dotted). 



Thus, we showed that the response of a droplet to com- 
pression is not harmonic, even to lowest order. It is not 
possible to obtain an expression for the potential as a 
function of £ in closed form, given the presence of elliptic 
integrals. It is possible, however, to test phenomenologi- 
cal expressions from our exact results. For this purpose, 
the divergence of the excess energy at h — > is not of par- 
ticular interest given the perspective of a droplet packing. 
We shall therefore concentrate on the droplet response to 
small compression. 

It is instructive to compare the forces predicted by the 
different estimates to the exact results. Figure || shows 
the ratio //£ for the exact solution, its lowest-order ex- 
pansion (Eq. |3l| ), the semi-circular surface model, and 
the truncated sphere model. A harmonic potential would 
be represented by a constant value in such a plot. For 
£ < 0.001, the exact solution sharply drops to zero as 
does its lowest order expansion. On the other hand, as 
£ goes to zero, the semi-circular surface model goes to 
(2 — 16/7T 2 ), showing that the potential becomes har- 
monic as £ — > 0. Finally, the truncated sphere model 
is found to be closer to a harmonic potential, having a 
smaller relative variation in //£. There are two essential 
characteristics of the exact force. The first one is the 
anharmonic behavior of the force at £ — > 0. This fea- 
ture is responsible for the logarithmic onset of G at ip c 
predicted for ordered emulsions || . The second essential 
characteristic is that //£ increases as the displacement 
increases. These important features, we believe, are re- 
sponsible for the characteristic positive curvature of G 
observed for disordered droplet packings. 
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FIG. 6. Logarithmic plot of the relative excess energy ver- 
sus the displacement function [(R/h) [i — 1]. A simple power 
law describes the data very well at low compression, as shown 
by the dotted line. The low-compression fits a slope of 2.1 
while the high-compression data joins the asymptotic limit of 
slope i. Curves are, top to bottom: truncated sphere ap- 
proximation (dot-dashed), semi-circular surface approxima- 
tion (dashed), and exact (solid). 



An approximate functional form for the droplet poten- 
tial would be useful in providing estimates of the elastic 
properties of compressed emulsions. For this purpose, 
the increase of //£ with the displacement can either be 
approximated by a power law £ a with a > 2 or by the 
addition of higher-order terms to a harmonic potential. 
The former is preferred however, both for its simplicity 
and because it vanishes as £ — > 0, similar to the exact 
solution. 

In view of investigating fitting functions, a logarith- 
mic plot of e in terms of [(R/h) 3 — 1] is presented in 
Fig. The numerical solution of Eq. |3l] can be performed 
with a precision larger than the one optimally possible 
when solving directly for e (Eq. |22|). Therefore, the low- 
compression values of e shown on Fij^ | were obtained 
using the lowest-order expansion (Eq. [52 ) . There are two 
distinct behaviors separated by a crossover section. By 
approximating the droplet by a compressed cylinder, it 
can easily be shown that e ~ 1/h as h — > 0. This asymp- 
totic limit is represented by the dotted line of slope 1/3 
at high abscissa values. At low compression, the data are 
relatively well described by a power law of the form £ Q or 
) 3 ~~ I]"" ^ e f° un d that the latter provides a better 
fit over a wider range of h. It also has the advantage of 
simplifying to the form (ip — <yS c ) Q for space-filling poly- 
hedra. Moreover, the two forms are equivalent for small 
£ since [(-f ) 3 — 1] ~ 3£. As we shall see in the next sec- 
tion, a fitting form will prove useful when characterizing 
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FIG. 7. The effective exponent characterizing the increase 
of the relative excess surface energy dlne/dlnx, where x 
is either the displacement £ (solid), or a function of it 
[1/(1 — — 1] (dashed). In both cases, the maximum value 
of £ is 2%. 

the behavior of the droplet for different configurations of 
compressing planes. 

As can be seen by comparing the curves of Fig. ^ at 
low compressions, a power law describes the data very 
well, down to machine precision values at e ~ 10 -8 . The 
fitting curve (dotted) has an exponent a = 2.1 and is ob- 
tained from fitting the data over a range covering more 
than an order of magnitude. The exact solution is con- 
trasted with the other estimates which show a harmonic 
behavior (a = 2) at small £. At larger £, the power 
law overestimates the excess energy but this occurs at 
£ sa 0.2. Nevertheless, the range of £ over which a power 
law describes the data very well is more than two decades. 

To investigate the validity of a power-law functional 
form for e at very small compressions, we simplify Eq. |3l] 
by taking one derivative with respect to £ and integrating 
/(£) to obtain 



£ = -^/M//4). 



(33) 



This allows us to define an effective exponent for very 
small compression 



dln£ = £/ = 21n(//4) 
dln£ e ln(//4) + l/2' 



(34) 



characterizing the increase of the excess energy. The ex- 
ponent is obtained in terms off which can be numeri- 
cally converted to £ using Eq. |3J. Figure shows the £ 
dependence of the effective exponent for the two different 
fitting forms we propose here. For both cases, points on 



the curve are bounded by £ < 2%, since these results are 
derived from the lowest-order expansion of the energy. 
In contrast to what is suggested by Fig. ||, the exponent 
is not constant for very small compression, but slowly 
decreases, going to 2 only at £ — * 0. This notwithstand- 
ing, a power law with a fixed exponent provides a good 
approximation for the droplet potential over the range 
of values of £ typically encountered in compressed emul- 



IV. NUMERICAL RESULTS 

Except for special cases, as shown in the previous sec- 
tion, the excess surface energy of a deformed droplet 
can only be determined numerically. Given the Sur- 
face Evolver software written by Brakke ||, it is rela- 
tively easy to calculate the shape of a single droplet with 
minimum surface area (energy) under the constraint of 
fixed droplet volume and as a function of confinement. 
We have obtained results for various confining polyhe- 
dral cells, some of which are Wigner-Seitz cells associ- 
ated with standard packing structures. The surprising 
result of this study is that for a considerable range of 
compressions, the energy is well described by a power 
law, where the exponent depends on the number of faces 
of the confining cell. 

The starting state is a spherical droplet which is com- 
posed of either 3074 or 12290 vertices. Most of the results 
presented here are for the latter case. The initial droplet 
is confined inside a polyhedral cell having between four 
and twenty faces. In particular, we used a rhombic dodec- 
ahedron (face-centered cubic), a truncated octahedron 
(body-centered cubic) and a simple cube (simple cubic), 
all of which are space-filling polyhedra. In addition, we 
compressed the droplet within a tetrahedron, an octa- 
hedron, a pentagon dodecahedron, and an icosahedron. 
Simpler structures were used as well, namely compres- 
sion between a pair of parallel plates to compare to the 
results obtained in Sec. Ill, and between two perpen- 



dicular pairs of parallel plates. Results for two parallel 
plates are shown as the open circles in Fig. ^. As ex- 
pected, they are in agreement with the analytic calcula- 
tions down to very low compression where the change in 
surface area is comparable to the numerical uncertainty 
in SE. In all cases, the compressions leave the center of 
mass unchanged so that a distance h from the center can 
be defined. 

Figure || shows the calculated dependence of the rela- 
tive excess energy (Eq. |24|) per facet, e/n, on £ (Eq. ||) 
for different types of confinement cells, n represents the 
number of compressing planes which is equivalent to the 
coordination number for space-filling cells. Using a fit- 
ting function of the form 



e/n = C 



(35) 
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0.018 




TABLE I. Values of the parameters C and a of the 
power- law fit to the relative excess energy (Eq. ^) . The vari- 
ous configurations are ordered according to their coordination 
number n. Errors are in the last digit. 



configuration 


n 


C(1CT 3 ) 


a 


pair of | plates 


2 


7 


2.1 


two _L pairs of || plates 


4 


12 


2.1 


tetrahedron 


4 


14 


2.1 


cube (s.c.) 


6 


21 


2.2 


octahedron (low-ip b.c.c.) 


8 


27 


2.3 


rhombic dodecahedron (f.c.c.) 


i2 


55 


2.5 


pentagonal dodecahedron 


12 


54 


2.5 


icosahedron 


20 


69 


2.6 



0.004 
0.002 
0.000 



0.08 



0.10 



0.00 0.02 0.04 0.06 

FIG. 8. Relative excess surface energy per facet, Eq. |24|, 
for a single droplet as a function of the plate displacement, 
£ = (R — h)/R. Curves are, from top to bottom, n = 12 
(f.c.c), 8 (low compression b.c.c), 6 (s.c), 4 (tetrahedron), 2 
(2 plates). All curves are best fit to power law Eq. |35[ Note 
that for clarity, not all cases of Table H are shown here. 



we find that both the coefficient, C, and the exponent, 
a, depend on n, as indicated in Table [I]. As mentioned 
before, Eq. |3| is equivalent to C'((p — <p c ) a when the 
droplet is confined by space-filling polyhedra. As in the 
case of two parallel plates, this form fits the data better 
over a wider range than £ a . The data selected for the 
fit lie in the £ interval 2-6%. This range was chosen to 
avoid errors associated with low values of £. On the other 
hand, the upper value fixes the range we are interested 
in. Indeed, since £ = 1 — (tp c /(p)3 , a face-centered-cubic 
structure (f.c.c, ip c = n\/2/6) for instance, has £ « 9.5% 
in the dry foam limit ip = 1. 

Our results explicitly show that the response of a 
droplet to compression is a non-local phenomenon: the 
response depends on the number of planes used to com- 
press the droplet. In addition, one might expect it to 
depend on the relative positions of these planes, but it 
will be shown below that this effect is not very impor- 
tant. It is worth noting that since a(n) > 2 for all cases, 
this response is weaker than a harmonic interaction: the 
response is equivalent to having the spring constant go to 
zero as the distortion vanishes. We also note that a power 
law with an exponent which increases as the number of 
compressing planes increases is unphysical for small com- 
pressions as the energy curves for different n would cross. 
Thus, although the fits shown in Fig. |^ appear convinc- 
ing, the functional form should be viewed uniquely as a 
convenient way to mimic the response of a droplet. 

In view of investigating the dependence of the relative 
positions of the compressing planes on the droplet, we ob- 
tained data for different configurations having the same 



n. By comparing, say, the results of the f.c.c. structure 
to those of the pentagon dodecahedron, we see that the 
influence of the configuration on the fitting parameters 
is marginal compared to that of n. It therefore seems 
reasonable to conclude that, provided that the different 
facets of the droplet are not too close, the potential is 
only a function of the coordination number n. This in- 
formation is valuable for modeling the excess energy of a 
droplet in a disordered emulsion jlj] . A second point of 
interest is the saturation of the exponent a above n sa 12. 
Although the number of neighbors nearly doubles from a 
rhombic dodecahedron to an icosahedron, the exponent 
barely changes. 

In a similar way, by comparing the excess energy ob- 
tained by compressing a droplet in an octahedron to 
that obtained for the truncated octahedron, one can de- 
termine the importance of the second neighbors in the 
base-centered-cubic (b.c.c.) lattice. We find that second 
neighbors do not play a role for <p < 0.90 (£ ~ 6.3%). 
For ip ~ 0.93, the excess energy differences between the 
two is about 5%. Thus over most of the range of interest, 
only the eight nearest neighbors are relevant. 



V. OPTIMAL PACKINGS FOR FOAM 

Our numerical results also allow us to compare the 
energy of the various packings as a function of volume 
fraction. For a wide range of volume fractions, the low- 
est energy state of monodisperse emulsions is believed 
to be a packing of droplets in an f.c.c. structure. We 
neglect temperature effects in the following argument. 
For ip > ip rn 0.545, i.e. the melting concentration 
of monodisperse hard spheres |22| , up to ip < ip c = 
7r\/2/6 ~ 0.7405, i.e. close packing concentration, the 
droplets remain spherical in an f.c.c. structure as they 
can pack without touching. For ip > ip c , the lowest en- 
ergy state remains f.c.c. until it changes to a new struc- 
ture at ip* < 1. That such a crossover must exist is clear 
from previous work on dry foams. 

In the limit (p « 1, an emulsion becomes a biliquid 
foam which is structurally analogous to a dry foam. The 



8 



0.12 




0.00 



0.65 0.70 0.75 0.80 0.85 0.90 0.95 1.00 
<P 

FIG. 9. The scaled excess energy density versus the vol- 
ume fraction as extracted from a droplet in a b.c.c. (□) and 
fee. (•) lattice. We find a crossover from fee. to b.c.c. 
at tp* = 0.932(1). tp c = ty^3/8 ~ 0.6802 for b.c.c. and 
iry/5/6 ~ 0.7405 for fee. 



lowest-energy state of the latter has long been the object 
of study. For nearly one hundred years, it was thought 
that the Kelvin |23| structure based on a b.c.c. pack- 
ing of identical orthic tetrakaidecadedra gave an optimal 
space filling of cells of the same volume. The mechan- 
ical equilibrium requirements of the foam structure are 
contained in Plateau's rules (l(J. To satisfy these rules, 
Kelvin showed that the faces of each cell are slightly 
distorted p3| |: the six quadrilateral surfaces remain pla- 
nar but curve their edges while the eight hexagonal sur- 
faces, in addition to sharing the same curved edges, be- 
come non-planar surfaces of zero mean curvature. The 
reduction in surface area from the undistorted Wigner- 
Seitz cell is approximately 0.16% @||[. By using SE, 
Weaire and Phelan |^6| recently showed the existence of 
a packing having an energy lower than the Kelvin struc- 
ture. This space-filling packing consists of six 14-sided 
polyhedra and two 12-sided polyhedra, all of equal vol- 
ume, packed in a simple cubic cell. This new structure, 
which is based on the cubic clathrate structure, has a sur- 
face energy approximately 0.3% less than that of Kelvin. 
For comparison, the f.c.c. structure has an energy 0.7% 
higher than the Kelvin structure |2(| . 

Since we have studied only a single droplet in a Wigner- 
Seitz cell, we cannot address the stability of the Kelvin 
structure versus that of Wcairc-Phclan for tp < 1 . How- 
ever, we can study the crossover from the f.c.c. packing 
to the b.c.c. Kelvin-like packing as tp increases. In Fig. ||, 
we plot the (scaled) energy density, U, defined by 



U = E/V = 3tpe(a/R) 



(36) 



as a function of tp for a droplet in b.c.c. and f.c.c. 
structures. Here, E is the excess surface energy E = 
<j(A — 4nR 2 ). As expected, the energy density of the 
f.c.c. structure is the lowest over a wide range of tp. But 
at tp* = 0.932(1), we observe a crossover to the b.c.c. 
structure. This result is in agreement with a similar pre- 
diction by Kraynik |27j] . Note that this crossover value is 
obtained strictly from surface energy considerations and 
may be superseded by other stability criteria. Moreover, 
our calculation neglects some degrees of freedom since 
it strictly involves compressing a droplet inside a polyhe- 
dral cell. This is more restrictive than compressing many 
droplets caged by their neighbors, and thus we are in fact 
obtaining an upper bound on the energy. However, since 
the energy of a drop in the shape of a truncated octahe- 
dron is only 0.16% higher than the orthic tetrakaidecad- 
edra of the Kelvin structure in the dry foam limit, the 
corrections to the crossover value resulting from the sur- 
face curvature of the films between the droplets should 
be small. Note that the crossover point is just above the 
volume fraction where second neighbors begin to play a 
role in the energy of a compressed droplet in the b.c.c. 
phase. Thus, while the six second neighbors contribute to 
increasing the deformation of a droplet and thereby the 
surface area, the overall energy of the b.c.c. structure is 
still lower than that of the f.c.c. structure for tp > 0.932. 
As we shall see in the next section, second neighbors in 
the b.c.c. structure are required for stability under shear 
deformation. 

The little difference in energy shown in Fig. ^| is not 
specific to b.c.c. and f.c.c. pairs. In fact, it is interesting 
to compare the energy per droplet obtained for the reg- 
ular pentagon dodecahedron (r.p.d.) to those obtained 
for the b.c.c. and f.c.c. lattices. Since the r.p.d. is not a 
space-filling object, tp represents the ratio of droplet to 
cell volume for this case. If one compares the energy per 
droplet as a function of the volume ratio, one finds that 
the results for the r.p.d. are almost indistinguishable from 
that of a truncated octahedron, i.e. the b.c.c. unit cell. 
In particular, results for the r.p.d. are lower than that 
for f.c.c. at high tp. This suggests that it may be favor- 
able for a compressed monodisperse emulsion to generate 
short-range icosahedral order such as it exists in metallic 
glasses. This implies that as an emulsion is compressed, 
it is likely to undergo a glass transition reminiscent of 
that occurring in metallic glasses. Indeed, it is found 
to be very difficult to obtain an ordered monodisperse 
emulsion experimentally, indicative of the deep energy 
minimum of an amorphous state. 



VI. RESPONSE TO SHEAR 

We now turn to the response of a compressed droplet 
to an impressed shear strain. Since we are considering 
a single droplet, the scope of our method is limited to 
ordered lattices. All the space-filling confinement cells 
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FIG. 10. Typical calculated excess energy density of a 
droplet compressed in an f.c.c. lattice as a function of shear 
strain. The strain is measured byA 2 + 2/A — 3~ 3e 2 so that 
the slopes are half the values of shear moduli |3lj . Curves are 
for, top to bottom, ip = 0.787, 0.763, 0.746. 

investigated here are related to the cubic point group 
and consequently, the associated elastic constant tensor 
is composed of three independent quantities. It is not 
possible for us to determine these constants individually. 

A useful method for applying a shear strain in pres- 
ence of periodic boundary conditions is to use an iso- 
choric uniaxial compressional-cxtcnsional strain. The ef- 
fective shear modulus can then be extracted from the 
resulting change in surface energy. The isochoric uniax- 
ial strain consists in an extension of the Wigner-Seitz cell 
of A = 1 + e in the, say, z direction, associated with a 
compression of A~2 in the perpendicular plane. In prin- 
ciple, different orientations of the Wigner-Seitz cell with 
respect to the reference z axis would lead to effective 
measures of the modulus resulting from different combi- 
nations of the elastic constants. However, considerations 
of the alignment of the facets through periodic boundary 
conditions put restrictions on the possible orientations of 
the strains. Or, said differently, some orientations of the 
uniaxial strain involve forces that are not normal to the 
faces of the Wigner-Seitz cell, and it results that some 
facets do not lie on the cell boundary 0. 

Here, the chosen orientation is along any of the three 
natural directions of a unit cell. For this particular 
choice, the strain tensor is written [§8| 

7mn = 2<W( A m ~ !)' ( 37 ) 

where S mn is the Kronecker delta, Ai = A2 = A~i, and 
A3 = A = 1 + e. e may be thought of as the magnitude 
of the small linear strain. The energy density U of the 



structure is given by 

where we use the summation convention. E is a fourth- 
order tensor which can be displayed as a 6 x 6 coefficient 
matrix using Voigt notation P^9|. When a small iso- 
choric uniaxial strain is applied, the energy density is 
found to vary as 

U=\Ge\ (39) 

where the effective shear modulus G is the combination 
of the elastic constants (Cu— Cvz)/2. Our definition of G 
reduces consistently to the corresponding Lame constant 
when considering an isotropic material. 

Figure [To] shows a typical result for the variation of 
surface area as a function of shear strain for an f.c.c. 
lattice at different volume fractions. The strain is mea- 
sured as A 2 + 2/ A — 3 which is equivalent to 3e 2 for small 
strains Ipjjl- The slope of each curve is G/2. 

In FigT|ll|, the shear modulus is shown for the s.c. and 
f.c.c. lattices as a function of volume fraction. These re- 
sults compare well with the ones obtained by Buzza and 
Cates II using Morse and Witten's potential Q . Similar 
to what they found for a s.c. lattice, the shear modulus 
of the f.c.c. structure shows a smooth rise (rather than 
a jump) at ip c . As we already mentioned, this behavior 
is different from that observed experimentally for disor- 
dered monodisperse emulsions. 

The isochoric uniaxial strain is peculiar in that with 
A = 2 3 , it changes the metric tensor in such way to trans- 
form a b.c.c. lattice into an f.c.c. lattice with the same 
density J30|. The net energy required for this transfor- 
mation will depend on the volume fraction as we have 
seen in the previous section when comparing the energy 
of both lattices. For ip > ip c and up to some ip' < ip*, the 
b.c.c. structure was found to be unstable (G < 0) to an 
applied uniaxial strain. The existence of ip' derives from 
the stability of the b.c.c. lattice over the f.c.c. lattice for 
ip > tp*. Second neighbors are responsible for stabilizing 
the structure against shear strain and it is not surprising 
to find that <p' > 0.90, the point where second neighbors 
start touching the droplet. Our best estimate for ip', i.e. 
the volume fraction above which a b.c.c. lattice has a 
positive shear modulus for an isochoric uniaxial strain, is 
0.903(5). 

VII. DISCUSSION AND CONCLUSION 

We have shown that the droplet response to small com- 
pressions in three dimensions is not harmonic and that 
it depends on the number of neighbors. For droplets 
compressed between two parallel plates, the surface pro- 
file can be solved analytically for small deformations and 
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FIG. 11. Shear modulus for a droplet in s.c. (□) and f.c.c. 
(A) Wigner-Seitz confinement cells as a function of compres- 
sion. Errors are of the order of 10%. ip c = t/6 ~ 0.5236 for 
s.c. and tt\/2/6 ~ 0.7405 for f.c.c. 



integrated numerically for arbitrary deformations. How- 
ever for other geometries, the deformation can only be 
determined numerically. Using the Surface Evolver soft- 
ware 0] , we have determined the change in surface area of 
droplets compressed in a variety of confining polyhedra. 
Our results strongly suggest that in three dimensions, 
though not in two, the response is non-local. The data 
fit very well with a non-harmonic scaling form in which 
the parameters of the fit depend almost exclusively on 
the coordination number. This feature has been used 
elsewhere O] to build a pairwise inter-droplet potential 
in order to investigate the elastic response of disordered 
monodisperse emulsions. 

The response to shear was also investigated for ordered 
lattices and it was found that similar to previous stud- 
ies [|J , the volume-fraction dependence of the shear mod- 
ulus shows a rather sharp increase at the onset of droplet 
touching, ip c . Although the anharmonicity of the poten- 
tial certainly plays a role in the linear onset of G for 
disordered emulsions, it was shown here that it is not suf- 
ficient. The experimental results currently available [^]]|] 
are for disordered emulsions only. Ordered lattices are 
easier to study theoretically but the effects of disorder 
seem to be too important to permit comparison with ex- 
periment. We have shown elsewhere |l4| that both an- 
harmonicity and disorder are required to reproduce the 
experimental results. 
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